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Abstract 

The swimming of a deformable planar slab in a viscous incompressible fluid is studied on the 
basis of the Navier-Stokes equations. A continuum of plane wave displacements, symmetric on both 
sides of the slab and characterized by a polarization angle, allows optimization of the swimming 
efficiency with respect to polarization. The mean swimming velocity and mean rate of dissipation 
are calculated to second order in the amplitude of the stroke. The optimum efficiency depends 
on the ratio of viscosity and mass density of the fluid. For high viscosity a stroke is found with 
significantly higher efficiency than Taylor’s solution for a waving sheet. For low viscosity the 
efficiency is optimal for a nearly irrotational flow pattern. 
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I. INTRODUCTION 


Historically the theory of swimming and flying has developed in two quite different di¬ 
rections [l|. On the one hand there is the work on swimming and flying of a plate at high 
Reynolds number started by von Karman and Burgers and continued by Lighthill j^,|3] 
and Wu 1-0- The work is based on Euler’s equations for irrotational flow of a fluid with 
vanishing viscosity. The discontinuity in velocity potential at the rear of the body is resolved 
by a vortex street. On the other hand there is the work on swimming at low Reynolds num¬ 
ber based on Stokes equations for a viscous incompressible fluid. This work was started by 
Taylor with his model of a swimming sheet, and continued for bodies of different shapes 
by Lighthill [l^, Purcell 11|, and Shapere and Wilczek among others. 

Taylor’s work on a swimming sheet was extended by Reynolds [l^, who considered the 
effect of fluid inertia on the swimming velocity. His work was corrected by Tuck [l^, who 
calculated the mean rate of dissipation as well. For transverse plane wave displacements of 
given wavenumber and frequency the fluid flow to first order in the amplitude is irrotational 
in the limit of small viscosity, but the relation to the work of Lighthill and Wu is not clear. In 
the latter the swimming speed is treated as a parameter, and the thrust is calculated instead. 
The concept of thrust plays no role in the work of Reynolds and Tuck. The swimming of a 
sheet in the limit of small viscosity has been studied by Childress 1^| by the technique of 
matched asymptotic expansions. 

In recent work on swimming and flying of an assembly of rigid spheres [T^ we have 
developed a mechanical model for which the swimming velocity and mean rate of dissipation 
can be calculated in the whole range of fluid viscosity and mass density. In the model the 
mean thrust, equal and opposite to the mean drag, vanishes identically in periodic swimming. 
Further investigation into different models in the full range of viscosity and inertia would 
be helpful. In the following we study small amplitude swimming of a slab, which continues 


earlier work [1^ on the resistive limit. We begin by rederiving Tuck’s expressions for a sheet 
using the method of Felderhof and Jones [l9|. Subsequently the same method is used to 
study the swimming of a slab. 

The main purpose of the present study is to investigate the effect of fluid inertia on 
the efficiency of swimming for a simple model. Tuck’s extension of Taylor’s work on the 
swimming of a sheet showed that in that case fluid inertia has a negative effect in the 
sense that the swimming velocity for a stroke of given power is always less than for purely 
resistive Stokes flow. The same is true for swimming by longitudinal distortions in the plane 
of the sheet, also investigated by Tuck 1^. On the other hand, for a slab one can consider 


swimming by potential flow, where viscosity plays no role and the swimming is achieved 
purely by inertia and the effect of added mass. As we show below, such potential swimming 
is just as efficient as that of Taylor’s sheet or the squirming sheet in the Stokes limit. 

After a discussion of Tuck’s work on transverse waves on a sheet we study hrst the 
swimming of a slab by irrotational flow. The flow pattern has a different symmetry from that 
considered by Lighthill and Wu, the hrst order how potential and pressure being continuous 
across the slab. The stroke consists of a circular motion of points on the surface of the slab, 
as in gravity waves, and of opposite circular polarization on both sides. There is no viscous 
boundary layer. It turns out that the swimming efficiency, dehned as the ratio of swimming 
speed and mean rate of dissipation, divided by viscosity, is independent of both viscosity 
and huid mass density. It has the same value as for Taylor’s sheet. 

Subsequently we study a squirming slab with displacements in the plane of the slab, so 
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that the displacement wave is longitudinal. It turns out that the swimming efficiency for 
this stroke is always larger than that of Tuck’s solution for the transverse wave. In both 
cases the efficiency tends to Taylor’s value in the limit of large viscosity and to zero in the 
limit of small viscosity. The latter behavior implies that these linearly polarized waves give 
a poor description of swimming and flying at small dimensionless viscosity C 

The wingtip path for birds is close to elliptic, or sometimes a hgure eight |2l|. To obtain a 
better description of swimming and flying at small viscosity it is therefore natural to consider 
strokes with a forward motion in the form of an elliptically polarized wave. We consider 
in particular a slab with symmetric displacements on both sides, given by an elliptically 
polarized plane wave, characterized by a polarization angle a. The displacement vector is 
given by Eqs. (2.12) and (6.1). For a = 0 the solution is identical to the squirming flow 
solution mentioned above. For a = vr/d it corresponds to the potential flow solution, and 
for a = 7r/2 it is identical to Tuck’s transverse solution apart from a sign factor for the flow 
in the lower half-space. 

The polarization angle a can be used to optimize the efficiency for hxed wavenumber and 
frequency. It turns out that in the resistive limit the efficiency is maximal for a = vr/S, with 
efficiency a factor \/2 larger than for Taylor’s solution, as we found earlier [l^. For smaller 
values of the dimensionless viscosity ( the optimum angle is in the range vr/S < a < vr/d 
and leads to a value of the optimum efficiency between one and \/2 times that of Taylor’s 
sheet. In the inertia-dominated limit, characterized by ^ 0, the optimum corresponds to 

a = 71 /4 and the potential flow solution. The optimum value of the efficiency in this limit 
is identical to that of Taylor’s sheet. In Figs, d and 5 we present a sketch of the nature of 
the motion at the surface of the slab for this case. The calculation suggests that in other 
geometry as well, the potential flow solution without viscous boundary layer is of paramount 
importance in the limit of low viscosity. 


II. SWIMMING SLAB 

We consider a planar slab immersed in a viscous incompressible fluid of shear viscosity 
rj and mass density p. The slab at rest is bounded by two planes at distance 2d. We use 
Cartesian coordinates x, y, z such that the upper xz plane is at y = d and the lower at 
y = —d. The rest shape of the slab is denoted as Sq. We shall consider a prescribed time- 
dependent shape S{t) leading to swimming motion in the x direction. The distortions are 
decomposed into S+{t) for the upper surface and *S'_(f) for the lower surface. The fluid is 
set in motion by the time-dependent distortions of the slab. The flow velocity v{r,t) and 
the pressure p{r,t) in the rest frame of the slab satisfy the Navier-Stokes equations 

d'v dTJ^ 

+ (i^ ■ v)u] = T]V^v - Vp - P— , V ■ u = 0 , (2.1) 

where U{t) is the instantaneous swimming velocity. We assume that the distortions are 
periodic in time with period T = 2t 7juj. They give rise to a mean translational velocity 
U = f/Ca; in the x direction. The slab can be idealized to be inhnite in the x and 2 ; 
directions. The flow velocity tends to —U{t) and the pressure tends to the ambient pressure 
Po as y tends to ± inhnity. We assume that the distortions of upper and lower surface have 
a symmetry such that they do not give rise to a rotational velocity. 

The surface displacement is dehned as the vector distance 

^ = s'-s (2.2) 
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of a point s' on the displaced surface S(t) from the point s on the slab So- We decompose s 
into s+ = (x, d, z) for the upper surface and s_ = (x, —d, z) for the lower surface. The fluid 
velocity v{r,t) is required to satisfy the no-slip boundary condition 


v{s + ^{s,t)) = 


dt 


(2.3) 


The displacement is decomposed into for the upper and lower surface. We consider 

displacements which do not depend on z and take the form = (^±a;(cc, f), t), 0). As 

a consequence the flow velocity v and pressure p do not depend on 2 ;, and the problem is 
effectively two-dimensional in the x, y coordinates. 

We construct an approximate perturbative solution to Eq. (2.1) with boundary condition 


(2.3) by formal expansion of the flow velocity and the pressure in powers of 


V = V 1 + V 2 +p = Pq+Pi+P2 + .... (2.4) 

Correspondingly the translational swimming velocity is expanded as 

U = Ui + U2 + .... (2.5) 

The flow velocity at the surface is formally expanded as 

?;(s + ^(s,f)) =v{s,t) + + .... (2.6) 

With the aid of this expansion the boundary condition may be applied at the undisplaced 
surface. 

To hrst order the fluid equations of motion reduce to the linearized Navier-Stokes equa¬ 
tions 

P-^ =d^^Vi-Vpi- V-?;i = 0. (2.7) 

at dt 

The hrst order boundary condition is 


r;i(s) 


dt 


s E Sq. 


( 2 . 8 ) 


To second order the time-averaged how V 2 , P 2 satishes the inhomogeneous Stokes equa¬ 
tions [19 

pV%2 - = -F 2 , (2.9) 


with force density 


F2 = ■ V)i>i, 


( 2 . 10 ) 


where the overline indicates averaging over a period T = 27r/a;. The force density may be 
written as the divergence of a Reynolds stress tensor. The Stokes equations Eq. (2.9) must 
be solved with the boundary condition 


«2(s) =-(^ ■ V)vi |g, (2.11) 

as follows from Eq. (2.3). We shall consider situations where the value on the right is 
the same for y = ±d, the mean how velocity V 2 = depends only on \y\, and tends 
to a constant —U 2 e.x as y tends to ±cx). The prefactor U 2 can be identihed as the second 
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order mean swimming velocity of the slab. The force density F 2 is of the form F 2 = 
F2x{,y)e^ + F 2 y{ii)ey and the mean pressure p 2 depends only on y. 

It will be convenient to use complex notation. Thus we write for the displacement vector 
of the upper and lower surface 

( 2 . 12 ) 

with complex amplitude and with the understanding that the real part of the expression 
is used to get the physical displacement. Correspondingly the hrst order velocity Ui(r,t) 
and pressure pi{r,t) take the form 

vi{r,t) = Pi{r,t) = (2.13) 


The boundary value in Eq. (2.11) may be expressed as 


^2(5) = -^Re (^± ■ V)ui|,^^ 


y=±d' 


(2.14) 


with values taken on either the upper or lower side of the slab. 

The power required equals the rate of dissipation of energy in the fluid. This can be 
calculated from the work done per unit area against the stress cr = p{Vv + (Vu)^) — pi. 


The mean rate of dissipation per unit area is to second order jl9 


D 2 =-Re vl ■ (Ti ■ (2.15) 

We have taken account of a factor 2 since there are equal contributions from the upper and 
lower half-space. 

We dehne the dimensionless efficiency as 


E 2 = (2.16) 

We used the prefactor 4 so that this equals unity for the case of a sheet in the resistive limit, 
as considered by Taylor j^. 


III. TUCK’S SOLUTION 


Tuck 1^ generalized Taylor’s solution for an undulating thin sheet to the case of a fluid 
with inertia. Here we derive his results from Eqs. (2.9)-(2.11) and cast them in a more 
elegant form. We use the displacement vector Eq. (2.12) with 


= kle, 


(3.1) 


with real amplitude A. In Fig. 1 we show a sketch of the slab with displacements given by 
Eq. (3.1). The hrst order how velocity and pressure are given by Eq. (2.13) with 


k 




(y) = v^^iy) = 


pliy) = 


s = \ — 


tup 

V 


B = 


s — k 


-H, 


(3,2) 
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where the upper (lower) sign refers to the upper (lower) half-space. In complex notation the 
force density F 2 in Eq. (2.10) can be expressed as 


Fj = P('"l ■ 


(3.3) 


The X-component has the value 
= IrA^pu^Re 


^ fu — 5 


*'\i,=F(s+s*)(y=F<i) 

(3.4) 


The corresponding solution of Eq. (2.9) takes the form V 2 = V 2 x{y)ex with component 

v^{y) = -U 2 


1 .9 2 i-i 

— — A^pu^Re — -, , 

2p |fc — sp -|- s 


£_ g=F(fc+s)(?;T<i) ^ g=F(fc+s*)(?/=F<i) _ 


k 


k + s* 


s s* 


^^:{s+s*){y^d) 


(3 


From Eq. (2.14) we hnd the boundary condition 

1 


V2x{^d) = 2^^^ R® 


Hence we hnd the swimming velocity 


1 


U 2 = -^A^ook Re/(C) 


with dimensionless viscosity 

and the complex function 


c = 


cop 


f{z) = - + ^z{i-z). 


The expression in Eq. (3.7) can be shown to be identical to that derived by Tuck [15 


Re/(C) = I^. F = + + forC>0. 


(3.6) 

(3.7) 

(3.8) 

(3.9) 

1 

(3.10) 


The variable F tends to unity at large C and diverges as 1/y/^ at small C- We note the 
identities 


y-i - c = lyc F+ 2 ^ f ' 


(3.11) 


Tuck calls 1/C the Reynolds number. In Fig. 2 we plot Re/(C) as a function of C- For large 
C, the function tends to unity in agreement with Taylor’s result [^. For C 0 the function 
tends to 1/2. 
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We also consider the mean rate of dissipation. From Eq. (2.15) we find 


D 2 = A‘^T]U}‘^k ( 1 + Im 


c 


(3.12) 


The viscous stress tensor does not contribute. It follows from Eq. (3.11) that the expression 
is identical to that derived by Tuck jl5| . 


D2 = A^riuj^k{l +F). 

It tends to Taylor’s result 2A‘^7]u‘^k as ( tends to inhnity, and it behaves as 

D 2 ^ A'^ as ?7 —)■ 0, 


(3.13) 


(3.14) 


independent of the wavenumber. In Fig. 3 we plot the efficiency E 2 , dehned by Eq. (2.16), 
as a function of (. The efficiency tends to unity for large ( and vanishes in the limit C —?• 0. 
It is of interest to consider also the kinetic energy of flow. We hnd in total per unit area 


1C2 = A-P^ 


2k 


1 + V? Iuia/T+C^ 


= A 


^puj 


2k 


2 / 1 

1 + wt; ■ 


2F 


(3.15) 


The expression shows that the mean virtual mass per unit area depends on the dimensionless 
viscosity. There is an interesting relation to the expression for U 2 in Eq. (3.7). 

It follows from Eq. (3.2) that for small ( the hrst order flow velocity consists of a potential 
flow and a thin boundary layer. In the limit C 0 only the potential flow contributes to 
the kinetic energy. In this limit the swimming velocity tends to U 2 = —\A^ujk, the mean 
rate of dissipation D 2 vanishes, the mean kinetic energy tends to /C 2 = puj‘^/{2k), and 
the efficiency vanishes, 7^2 = 0. The mean kinetic energy is proportional to the mean thrust 
calculated in inviscid irrotational flow theory 13,0. In the present theory the concept of 
thrust does not appear. In the theory of inviscid flow end effects dominate and the swimming 
velocity is left undetermined. There are conceptual differences here that do not allow easy 
comparison. 


IV. POTENTIAL FLOW 

In this section we compare the above results for the Taylor sheet with a second solution 
where the slab undulates in such a way that the flow is irrotational. In Taylor’s solution 
the sheet undulates such that to hrst order the sheet is not extended. For the potential how 
solution considered below the upper and lower surface of the slab are extended to hrst order. 
We use the displacement vector Eq. (2.12) with complex amplitude vector 

^1 = ^i(^x±iey), (4.1) 
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corresponding to a circularly polarized wave. In Fig. 4 we show the displacements for a 
period. The hrst order flow velocity and pressure are given by Eq. (2.13) with 


THvTd) 


vl±{y) = 


v‘Jy) = ±uT e 


THvTd) 


(4.2) 


Note that the hrst order pressure takes the same value on both sides of the slab, unlike Eq. 
(3.2). From Eq. (2.14) we hnd the boundary condition 

v^{±d) = k. (4.3) 

Hence we hnd the swimming velocity 

t/2 = ^A^uk, (4.4) 


independent of viscosity. From Eq. (2.15) we hnd for the mean rate of dissipation 

D 2 = 2A‘^Tju‘^k. (4.5) 

The pressure does not contribute. For the efficiency dehned in Eq. (2.16) we hnd E 2 = 1, 
independent of viscosity. For the kinetic energy of how per unit area we hnd 


IC 2 


A 


2 ^ 

4k ’ 


(4.6) 


again independent of viscosity. 

The expressions derived above show that the potential how solution is quite simple, as 
it is for a distorting sphere [ 2 ^. The swimming is more efficient than for Tuck’s solution 
of Sec. Ill for any value of the dimensionless viscosity (, as shown in Fig. 3. There is no 
viscous boundary layer in the hrst order how. The mean swimming velocity is the same as 
for Taylor’s sheet at the same mean rate of dissipation, but the mean kinetic energy of how 
is one third of that of Taylor’s sheet. 


V. SQUIRMING SLAB 


We consider also a squirming slab with equal plane wave displacements in the upper and 
lower plane of the form Eq. (2.12) with complex amplitude vector 


= kle^- 


(5.1) 


In the limit of high viscosity the mean swimming velocity and the mean rate of dissipation 
are the same as for Taylor’s sheet for the same amplitude of displacement A, but with 
swimming velocity in the opposite direction l^.jj^. 

The hrst order how velocity and pressure are given by Eq. (2.13) with 


-'xit 


(y) = - I n^±(i/) = ±Buj 

k 


vUv) = 


S = 4/— 


lup 

5 

7 ] 


B = 


k — s 


-A. 


(5.2) 
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The x-component of the mean force density in Eq. (3.3) has the value 


1 .9 9^ ‘ik 

F,^ = -A Re 


(/t+s)seTl'=+*)(»T‘‘4«:(«:+s*)eT<'=+‘9(»T<i)_s(s+s")eT(>+»*)(»T<i)V 

(5.3) 


The corresponding solution of Eq. (2.9) takes the form V 2 = V 2 x{y)^x with component 
v^{y) = -U 2 


ts 


— —A^pca^Re ,, , 

2p |fc — sp + s 


^zf{k+s){y^d) 


^_ g=F(fc+«*)(?/=F<i) _ ^ g=F(^i+s*)(?/=F<i) 


k + S 


S + S* 


From Eq. (2.14) we hnd the boundary condition 


V2x{,±d) = ~A^uk. 


Hence we hnd the swimming velocity 


U2 = (I - VC imvrrcf 


1 .2 , 3F- 1 

= -A^uk -. 

2 2F 


The second form was given by Tuck [15 


From Eq. (2.15) we hnd for the mean rate of dissipation 

= A^r]u^k(^l + -^Re^/-i + = A^r]u^k{l + F). 


The mean kinetic energy of how is 


ic2 = a 




ikF 


(5 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


In Fig. 3 we plot the efficiency E 2 , dehned by Eq. (2.16), as a function of (. The efficiency 
tends to unity for large ( and vanishes in the limit C 0 . 


VI. EFFICIENT SWIMMING 

It is evident from Fig. 3 that of the three modes considered so far, the potential how 
solution of Sec. IV leads to the most efficient swimming. In this section we show that for 
any value of the dimensionless viscosity ( = yk"^/ {oop) one can hnd a mode which is even 
more efficient. We consider intermediate modes depending on a polarization angle a which 
can be varied in such a way that the efficiency E 2 is optimized. The modes considered so 
far correspond to particular values of a. 

We consider plane wave displacements of the upper and lower plane of the form Eq. (2.12) 
with complex amplitude vector 

= Hcosa 63 . ± iHsina Cy, ( 6 . 1 ) 
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corresponding to an elliptically polarized plane wave. For a = 0 this corresponds to the 
squirming displacement of Sec. V. For a = 7r/4 it corresponds to the potential flow solution 
of Sec. IV. For a = ti/2 it corresponds to Tuck’s solution of Sec. Ill with displacement 
and flow in the lower half-space differing by a minus sign. The sign change corresponds to 
a shift of the flow pattern by half a wavelength in the x direction and does not affect the 
mean second order force density or flow velocity. The first order flow velocity and pressure 
are found by linear combination of the expressions in Eqs. (3.2) and (5.2). 

The x-component of the mean force density in Eq. (3.3) takes the form 

F 2 x = - sin2a) - 1 -/2x2 cos2a], (6.2) 

with coefficients given by 




i 

2\k-s\^ 

i 

2\k — sP 


(^s\k + s|2e=F(^+dfo=F<^) + k\k + ^| 2 gT(fc+s*)( 2 /T<i) _ 2ks{s + , 

(s{k + s){k- + k{k -s){k + . (6.3) 


Evidently for a = 7r/4 the expression in Eq. (6.2) vanishes, in accordance with Sec. IV. 
From Eq. (6.2) one evaluates the second order flow component v^{y) as in Eqs. (3.5) and 
(5.4). From Eq. (2.14) we hnd the boundary condition 


V 2 x{Fd) = --A^uj Re [/c — s -|- (/c -|- s) sin 2a + (k + s) cos 2a]. 


Hence we hnd the swimming velocity 


[/, = -A^uk 


1 1 {1 1 , 

-— -|- —I -— sin 2a -|- cos 2a 

2 2F V 2 2F 


Similarly we hnd for the mean rate of dissipation 

D 2 = A^rju'^k [1 -I- F -1- (1 — F) sin 2a], 
and for the mean kinetic energy of how 


/Co = H 


2k 


sin^ a -I-— (1 — sin 2a) 

2F 


(6.4) 


(6.5) 


( 6 . 6 ) 


(6.7) 


The mean rate of dissipation D 2 has a minimum at a = 7r/4, corresponding to potential 
hrst order how, independent of the value of C,. The ratio IF 2 I/F 2 yields the efficiency F 2 
according to Eq. (2.16). For hxed value of C, this can be maximized by variation of the 
polarization angle a. In the limit of large C, the maximum is at a = tt/S and then the 
efficiency takes the value \/2- For large the mean kinetic energy of how /C 2 is minimal 
at a = tt/S. For decreasing C, the maximum of F 2 shifts to larger values of a, tending to 
7r/4 as (C —)■ 0. At all values of C, the maximum is larger than unity, and it tends to unity 
as a —)■ 7r/4. For optimal swimming at small C, the how held is nearly irrotational. In Fig. 
4 we plot the distortions of the slab for purely irrotational swimming. We have assumed a 
symmetry between the upper and lower half space, but the how pattern in the lower half¬ 
space can be shifted horizontally without change in swimming speed or rate of dissipation. 
In particular, the distortions of the slab look like those shown in Fig. 5 for a shift by half a 
wavelength. 
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VII. DISCUSSION 


The expressions derived in Sec. VI show that for symmetric plane wave displacements of 
a slab the swimming performance depends in a complex way on viscosity and mass density 
of the fluid and on the polarization angle characterizing the stroke. The dependence allows 
optimization of the swimming efficiency by variation of the polarization angle. In the resistive 
limit there is a stroke whose efficiency is larger by a factor \/2 than that of Taylor’s sheet. 
In the inertial limit the optimal solution is found for nearly irrotational flow. 

The calculation of mean swimming velocity and mean rate of dissipation has been per¬ 
formed to second order in the amplitude of displacements. Clearly it would be of interest to 
carry the calculation to higher order, but it may be expected that the present calculation 
provides a good picture of the qualitative trend. 

The calculation is purely kinematic in the sense that the surface displacements are pre¬ 
scribed. A more elaborate calculation would take account of the elastic properties of the 
slab. The swimming as a consequence of an actuating force density could then be studied in 
a dynamic context. Also it would be of interest to extend the calculation to a compressible 
fluid. 
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Figure captions 


Fig. 1 

Plot of the slab with plane wave distortion corresponding to Tuck’s solution discussed 
in Sec. III. A cross-section in the xy plane with a length of one wavelength is shown. The 
motion does not depend on the coordinate. 

Fig. 2 

Plot of the function Re/(C) given by Eq. (3.10). The mean swimming velocity in Tuck’s 
solution is proportional to this quantity, as shown in Eq. (3.7). 


Fig. 3 

Plot of the efficiency E 2 as a function of dimensionless viscosity C, = / (cap) for Tuck’s 
solution (solid curve), for the squirming slab (long dashes), and for the irrotational flow 
solution (short dashes). 


Fig. 4 

Plot of the slab with symmetric plane wave distortion corresponding to the potential flow 
solution discussed in Sec. IV. A cross-section in the xy plane with a length of one wavelength 
is shown. The motion does not depend on the coordinate. 


Fig. 5 

Plot of the slab with plane wave distortion corresponding to the potential flow solution 
discussed in Sec. IV with the wavepattern for the lower half-space shifted in the x direction 
by half a wavelength in comparison with Fig. 4. 
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FIG. 5: 
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